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Abstract 

In light of the latest neutrino oscillation data, we revisit the minimal scenario of type-I seesaw model, 
in which only two heavy right-handed Majorana neutrinos are introduced to account for both tiny 
neutrino masses and the baryon number asymmetry in our Universe. In this framework, we carry out a 
systematic study of the Frampton-Glashow-Yanagida ansatz by taking into account the renormalization- 
group running of neutrino mixing parameters and the flavor effects in leptogenesis. We demonstrate 
that the normal neutrino mass ordering is disfavored even in the minimal supersymmetric standard 
model with a large value of tan/?, for which the running effects could be significant. Furthermore, it is 
pointed out that the original scenario with a hierarchical mass spectrum of heavy Majorana neutrinos 
contradicts with the upper bound derived from a naturalness criterion, and the resonant mechanism 
with nearly-degenerate heavy Majorana neutrinos can be a possible way out. 
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1 Introduction 


Neutrino oscillation experiments in the past two decades have revealed that neutrinos are actually massive 
particles and lepton flavors are significantly mixed [1], In order to accommodate tiny neutrino masses, 
one can go beyond the minimal Standard Model (SM) and introduce three right-handed neutrinos IV iR 
(for i = 1,2,3), which are singlets under the SU(2) L x U(1) Y gauge group of the SM. The most general 
gauge-invariant Lagrangian relevant for lepton masses and flavor mixing can be written as 

- An = hXl HE i l + kXuHN K + l -NlM K N K + h.c. , (1) 

where and H = ia 2 H* denote the left-handed lepton and Higgs doublets, respectively, while E R the 
right-handed charged-lepton singlets. In addition, Y) and Y v stand respectively for the Yukawa coupling 
matrices of charged leptons and neutrinos, and M R is the Majorana mass matrix for right-handed neutrino 
singlets. After the Higgs field acquires its vacuum expectation value ( H) = v ~ 174 GeV and the gauge 
symmetry is spontaneously broken down, the charged-lepton mass matrix is given by M , = Yju, while 
the Dirac neutrino mass matrix is M D = Y v v. Since the Majorana mass term for right-handed neutrino 
singlets is not subject to the electroweak gauge symmetry breaking, the absolute scale of M R could be 
much higher than the electroweak energy scale A EW ~ 100 GeV. Therefore, in the low-energy effective 
theory with heavy Majorana neutrinos integrated out, the mass matrix of three light neutrinos is given by 
the famous seesaw formula M u ~ Given 0(M- D ) ~ A E w, one can obtain neutrino masses 

at the sub-eV level if 0(M R ) ~ 10 14 GeV is close to the scale of grand unified theories A GUT ~ 10 16 GeV. 
In this canonical seesaw model El El n El Ej, the lightness of ordinary neutrinos can be ascribed to the 
heaviness of right-handed Majorana neutrinos. Moreover, the mismatch between the diagonalization of 
M l and M v leads to lepton flavor mixing. 

In the basis where both the charged-lepton mass matrix = diag {m e ,m fi ,m T } and the mass matrix 
of heavy Majorana neutrinos M R = diag {M 1 , M 2 , M 3 } = M R are diagonal, the neutrino mass spectrum 
and lepton flavor mixing are determined by the effective neutrino mass matrix M u = —M D M R 4 M R , which 
can be diagonalized as M u = U ■ diag {m 1 ,m 2 ,m 3 } ■ t/ T with U being the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrix (7» StEj- Therefore, in order to obtain any predictions for the low-energy 
observables, one has to know the flavor structure of M D , which is completely unconstrained in the generic 
seesaw model. Generally speaking, there are two different guiding principles towards seeking a solution to 
this problem, namely, flavor symmetry and minimality: 

• In the first approach, discrete or continuous flavor symmetries are imposed on the generic Lagrangian 
in Eq. m, and all the SM fields are assigned into proper representations of the symmetry groups. Due 
to the required symmetries, the Yukawa coupling matrices are not arbitrary any more. It has been 
demonstrated that discrete flavor symmetries can be implemented to successfully predict interesting 
lepton flavor mixing patterns, which are well compatible with the latest neutrino oscillation data. 
For recent reviews on this topic, see Refs. mmm- Although this scenario is very attractive in 
the first place, it actually suffers from the involvement of many new scalar fields that are needed in 
order to achieve the desired flavor structures of Yukawa coupling matrices. As a consequence, it is 
generally difficult to verify or disprove a flavor-symmetry model experimentally. 

• In the second approach, the number of model parameters is intentionally reduced to a level, beyond 
which the model would immediately run into contradictions with current experimental observations. 
The minimality of a model, in the sense of minimal number of free parameters, can be regarded as 
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an Occam’s razor Banal. One practical way of reducing free parameters is to simply take some 
Yukawa matrix elements to be zero. The physical essence of texture zeros actually reflects that some 
elements in a Yukawa coupling matrix are highly suppressed when compared to the other elements, 
or they are irrelevant to fermion mass spectra and flavor mixing. For instance, the texture zeros 
turn out to be very useful to establish a relationship between small flavor mixing angles and strong 
mass hierarchy in the quark sector Banana- As shown by Weinberg in Ref. na, the texture 
zeros in two-generation quark mass matrices lead to a successful prediction for the Cabbibo angle 
9 C = \Jm d /m s ~ 0.22, where the running mass of down quark m d = 2.82 MeV and strange quark 
m s = 57 MeV are evaluated at M z = 91.2 GeV [18, 19]. In the same spirit, more than ten years 
ago, Frampton, Glashow and Yanagida proposed a minimal scenario of seesaw models, in which only 
two right-handed neutrinos are introduced and two elements of the Dirac neutrino mass matrix M D 
are assumed to be vanishing [20]. In this case, M D becomes a 3 x 2 matrix, and can be explicitly 
written as 


M n = 



(2) 


V 


where a, 6, a' and b' are in general complex. There are totally fifteen possible patterns of M D with 
two texture zeros in different positions, and we shall examine all of them in the following section. The 
number of texture zeros in M D cannot be further increased, otherwise the model will be in conflict 
with three nonzero flavor mixing angles, as measured in neutrino oscillation experiments |2a EH 1221 
[23] . On the other hand, the seesaw model with just one heavy right-handed neutrino does not work, 
since there will be two massless ordinary neutrinos that have already been excluded. Hence, the 
scenario of two heavy right-handed neutrinos together with the Frampton-Glashow-Yanagida (FGY) 
ansatz like that in Eq. (J2J) is the minimal version of type-I seesaw model, which will be called the FGY 
model hereafter. One can immediately verify that neutrino mass spectrum and leptonic CP-violating 
phases are calculable from the observed three neutrino mixing angles and two neutrino mass-squared 
differences [23], implying a complete testability of the model in future neutrino experiments. It is 
worthwhile to stress that this minimal scenario emerges when one right-handed Majornana neutrino 
is much heavier than the other two and decouples from the theory, or its Yukawa couplings to lepton 
and Higgs doublets are vanishingly small m- 


Another salient feature of the canonical seesaw model is to account for the baryon number asymmetry in 
our Universe via the leptogenesis mechanism [23]. In the early Universe, the temperature is high enough 
to thermally produce heavy Majorana neutrinos As the Universe cools down, the out-of-equilibrium 
and CP-violating decays of N t generate lepton number asymmetries, which will further be converted into 
the baryon asymmetry via nonperturbative sphaleron processes [25l 126] , Excellent reviews on leptogenesis 
can be found in Refs. |271 [281129]. 

In light of recent progress in neutrino oscillation experiments, we reconsider the FGY model and carry 
out a complete study with a focus on the currently unresolved problems, such as neutrino mass ordering, 
leptonic CP violation and the Majorana nature of neutrinos. The main motivation for such an investigation 
is two-fold. First, due to a minimal set of free parameters, the FGY model is quite predictive, so it is 
interesting to confront it with the latest global-fit results of neutrino oscillation data. A similar analysis has 
actually been done in Ref. m- Different from that work, we take into account the renormalization-group 
(RG) running effects of lepton flavor mixing parameters from the seesaw scale A gg , usually characterized 
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by the lightest heavy Majorana neutrino mass M,, to the electroweak scale A EW . Second, in the previous 
work, a strong mass hierarchy M 2 is always assumed, and a narrow range of heavy neutrino masses 

M 1 ~ 5 x 10 13 GeV is derived by requiring a successful leptogenesis mechanism to explain the cosmological 
matter-antimatter asymmetry. But such a large mass scale in the theory causes the naturalness or fine- 
tuning problem on the one hand [3Ql EU |32l |33l [34] , and the gravitino overproduction problem if the model 
is supersynmretrized on the other hand [35]. Therefore, we are motivated to go beyond the hierarchical 
limit, and consider both mild mass hierarchy and a nearly-degenerate mass spectrum of heavy Majorana 
neutrinos. Only with careful studies of RG running effects and general mass spectra of heavy Majorana 
neutrinos can we really test the FGY model. 

The remaining part of our paper is organized as follows. In Section 2, phenomenological implications 
of the FGY model are explored and confronted with current neutrino oscillation data. We also consider 
the RG running effects of neutrino mixing parameters, and specify the allowed regions of the parameter 
space at the low-energy scale. Only four out of fifteen patterns of the Dirac neutrino Yukawa coupling 
matrices are found to be compatible with neutrino oscillation data, and only the inverted neutrino mass 
ordering is allowed. Section 3 is devoted to the generation of baryon number asymmetry via leptogenesis, 
where we also discuss the impact of lepton flavor effects and non-hierarchical mass spectrum of heavy 
Majorana neutrinos. The flavor structure of four viable patterns leads to a non-vanishing CP asymmetry 
in one specific lepton flavor. We point out that a nearly-degenerate mass spectrum of heavy Majorana 
neutrinos is required to explain the baryon number asymmetry, and simultaneously avoid huge radiative 
corrections to the light Higgs boson mass. Finally, we summarize our main conclusions in Section 4. 

2 Neutrino Masses and Flavor Mixing 

We start with neutrino mass spectrum and flavor mixing parameters in the type-I seesaw model with only 
two right-handed heavy Majorana neutrinos. After some general remarks, we proceed to introduce the 
FGY ansatz and explore its phenomenological implications. The RG evolution of neutrino masses and 
mixing parameters is considered when we confront the FGY ansatz with low-energy neutrino oscillation 
data. Finally, the model parameters relevant for leptogenesis at the high-energy scale are determined. 

2.1 General Remarks 

In the basis where both the charged-lepton mass matrix M l and the heavy Majorana neutrino mass 
matrix M R are diagonal, the diagonalization of the light neutrino mass matrix M v = via 

M v = UM^U 1 - gives us neutrino mass eigenvalues M u = diag{m 1 ,m 2 ,m 3 } and the PMNS matrix U. 
Since only two right-handed neutrinos are introduced and their mass matrix M R is of rank two, it is 
straightforward to verify that the rank of effective neutrino mass matrix M u is two. As a consequence, 
the lightest neutrino must be massless. In the case of normal mass ordering (NO) with m 1 = 0, we get 
m 2 = \J Am 21 and m 3 = y A m 3l . In the case of inverted mass ordering (10) with m 3 = 0, we have 
m l = \/|Am 3 2 | — Am 21 and m 2 = y| Am^l- The neutrino mass-squared differences Am 21 = m 2 — m\ 
and A m 31 = m 3 — (or Am| 2 = m 3 — m\) can be measured in neutrino oscillation experiments in the 
case of NO (or 10). At present, however, it is unclear whether neutrino mass ordering is NO or 10. The 
ongoing long-baseline accelerator experiments T2K [36] and NOz^A mi the forthcoming medium-baseline 
reactor experiments JUNO [38] and RENO-50 [39], and the future huge atmospheric neutrino experiment 
PINGU [50] will provide a definitive answer to this question. 
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Table 1: The best-fit values, together with the la, 2a and 3<r intervals, for three neutrino mixing angles 
{^ 12 )^ 13 )^ 23 }) t wo mass-squared differences or Amj 2 } and the Dirac CP-violating phase 5 

from a global analysis of current experimental data m • Two independent global-fit analyses can be found 
in Refs. mms], which are in perfect agreement with the results presented here at the 3 a level. 


Parameter 

Best fit 

lcr range 

2 a range 

3d range 

Normal neutrino mass ordering (m 1 < m 2 < m 3 ) 

$ 12 /° 

33.48 

32.73 — 34.26 

31.98 — 35.04 

31.29 — 35.91 

CO 

O 

8.50 

8.29 — 8.70 

8.08 — 8.90 

7.85 — 9.10 

O 23/ 0 

42.3 

40.7 — 45.3 

39.1 — 48.3 

38.2 — 53.3 

5/° 

306 

236 — 345 

0 — 24 © 166 — 360 

0 — 360 

Am^/flO" 5 eV 2 ] 

7.50 

7.33 — 7.69 

7.16 — 7.88 

7.02 — 8.09 

Am^/flO” 3 eV 2 ] 

+2.457 

+2.410-1-2.504 

+2.363 — +2.551 

+2.317-1-2.607 


Inverted neutrino mass ordering (m 3 < m 1 < m 2 ) 


$ 12 /° 

33.48 

32.73 — 34.26 

31.98 — 35.04 

31.29 — 35.91 

CO 

0 

8.51 

8.30 — 8.71 

8.09 — 8.91 

7.87 — 9.11 

023/ 0 

49.5 

47.3 — 51.0 

45.1 — 52.5 

38.6 — 53.3 

5/° 

254 

192 — 317 

0 — 20 0 130 — 360 

0 — 360 

Au4/[10- 5 eV 2 ] 

7.50 

7.33 — 7.69 

7.16 — 7.88 

7.02 — 8.09 

Am| 2 /[10” 3 eV 2 ] 

-2.449 

-2.496 -2.401 

-2.543 -2.355 

-2.590 -2.307 


Furthermore, the PMNS matrix in this minimal model can be parametrized via three mixing angles 
{<? 12 , $ 13 , 023 }) one Dirac-type CP-violating phase 5 and one Majorana-type CP-violating phase a, namely 


U = 

^ c 12 c 13 s 12 c 13 s 13 e 1<5 ^ 

— s 12 c 23 — c 12 s 13 s 23 el15 + c 12 c 23 — s 12 s 13 s 23 el<5 C 13 S 23 

| 

(l 0 o\ 

0 e ia 0 


V+ s 12 s 23 — c 12 s 13 c 23 el<5 ~ c 12 s 23 — s 12 s 13 C 23 el<5 C 13 C 23 ) 


T— I 

O 

O 


(3) 


where c- = cos $ - and -s- = sin 9-- have been defined for ij = 12,13,23. While three mixing angles have 
been determined with reasonably good precision from oscillation experiments, there is still no significant 
evidence for a nontrivial Dirac CP-violating phase. In Table 1, the latest global-fit analysis of neutrino 
oscillation parameters has been presented. One can observe that the best-fit value of Dirac CP-violating 
phase is 5 = 306° for NO and <5 = 254° for 10, but it becomes arbitrary at the 3 a level. The proposed 
neutrino super-beam experiments and neutrino factories are able to probe 5 down to a few degrees [44| . 

Since there is one massless neutrino, we have only one Majorana CP-violating phase a. The observation 
of neutrinoless double-beta decays is the unique and feasible way to establish that neutrinos are Majorana 
particles, i.e., they are their own antiparticles [35]. The decay rate depends on the effective neutrino mass 
defined as mop = \U^ 1 m 1 + t/ 2 2 m 2 + D 2 3 m 3 1, where U ei for z = 1,2,3 denote the elements in the first row 
of the PMNS matrix U. More explicitly, 

\J Am 31 cos 2 $ 13 [£ 2 sin 4 $ 12 + tan 4 $ 13 + 2£ sin 2 $ 12 tan 2 $ 13 cos 2 (a + 5)] 1 for NO ; 


m 


■pp 


= 


\J\ Am 32 | cos 2 $ 13 cos 2 $ 12 [C 2 + tan 4 $ 12 + 2( tan 2 $ 12 cos 2<r] 


1/2 


(4) 

for 10 , 


where £ = m 2 /m 3 and £ = rn x /m 2 . Now that neutrino masses are completely fixed by two mass- 
squared differences, we can get £ = \JAm^i/\J Am 31 ~ 0.175 and ( = \J 1 — Am 21 /|Am 32 | ~ 0.985 
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by using the best-fit values of neutrino mass-squared differences in Table 1. Notice that the relation 
£ 2 ~ 1 — (j 2 ~ \/2sin 2 9 13 ~ 0.03 holds as an excellent approximation. The exact value of Tripp depends on 
the Majorana CP-violating phase a in the 10 case, and a combination of two unknown CP-violating phases 
a and <5 in the NO case. However, it is straightforward to find out the lower and upper limits |.4Bi l3ZU3S] • 
For NO, we get 

yj cos 2 9 l3 (£ sin 2 9 12 — tan 2 9 13 ) < mpp < \J A m 3l cos 2 9 13 (£ sin 2 9 12 + tan 2 9 13 ) , (5) 

leading to mpp G [1.5, 3.7] meV with the help of the best-fit values in Table 1. For 10, we arrive at 

\J\/S.m\ 2 \ cos 2 9 l3 cos 2 0 l2 (£ — tan 2 9 l2 ) < rripp < yJ\Am% 2 \ cos 2 9 l3 cos 2 9 12 (C + tan 2 9 12 ) , (6) 

implying mpp € [18,48] meV with the best-ht values as inputs. As the future neutrinoless double-beta 
decay experiments are able to reach a sensitivity of about 20 meV [45], the IO case seems to be more 
encouraging and phenomenologically interesting. Moreover, in this minimal seesaw model, the observation 
of neutrinoless double-beta decays may also pin down the unique Majorana CP-violating phase a via 
Eq. ((4]), as long as the other mixing parameters can be well measured in neutrino oscillation experiments. 


2.2 The Frampton-Glashow-Yanagida Ansatz 


Although neutrino mass spectrum can be fixed by the observed neutrino mass-squared differences in the 
minimal seesaw model, three mixing angles and two CP-violating phases are in general arbitrary. Further 
restrictions on the flavor structure can induce testable correlations among low-energy observables. In the 
full theory above the seesaw scale A gs , relevant parameters are the Dirac neutrino Yukawa coupling matrix 
Y u and heavy Majorana neutrino masses {M 1 ,M 2 }. If two elements of Y u are vanishing [20], there are 
fifteen logically possible patterns, which can be categorized into three classes: 


• Case A Two texture zeros are located in the same row, namely, ( Y v ) ai = (Y^)^ = 0 with i / j. 
There are only three patterns: 


( 0 0^ 


/x x\ 


/x x\ 

X X 

, A 2 : 

0 0 

, A 3 : 

X X 

\x xj 


U x) 


1° V 


(7) 


where the cross ‘x’ denotes a nonzero matrix element. 


• Case B - Two texture zeros are located in different columns and rows, namely, (Y u ) ai 
with a / /3 and i / j. There are six patterns: 


/° x\ 


7 o x ^ 


(x x\ 

x 0 

, B 2 : 

X X 


0 x 

\x x) 


lx o) 


l x 0 / 


B 4 : 

7x 0^ 

0 x 

, B 5 : 

7x 0^ 

X X 


X X 

o x^ 


lx x) 


1° xj 


xJ 


where the patterns B 4 5 6 are derived from B x 2 3 by exchanging two columns. 


(Y u ) 0j = 0 


( 8 ) 
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Case C - Two texture zeros are located in the same column, namely, (Y u ) ai = (Y^C- = 0 with 
a 7 ^ (3. There are six patterns: 



/° x\ 

( 0 x ^ 


/x x\ 

Ch 

0 x 1 j C 2 ■ 

X X 

. C 3 ‘ 

0 x 


lx x/ 

1° xj 


1° x/ 


(9) 



fx 0\ 


fx 0\ 


/x x\ 

C 4 : 

0 ^ 

XX 

. C 5 : 

X X 

0 X 

. c 6 : 

x 0 

lx 0 / 


where the patterns C 4 5 6 can be obtained from C 4 2 3 by exchanging two columns. 

It is worth pointing out that the patterns in each class can be related by the elementary transformations, 
i.e., the 3x3 elementary matrices (for ij = 12, 23,13) and the 2x2 elementary matrix Q. The action 
of Vjj from left (or right) induces an exchange between z-th and j-th rows (or columns), and likewise for 
Q. With the help of and Q, one can change the positions of texture zeros. For instance, we have 
Y u ( A 2 ) = 7 : ’ 12 Y' jy (A 1 ) and Y u ( A 3 ) = V 13 Y U (A 1 ). In a similar way, one can prove that all the patterns in 
Case B can be obtained from Y u ( B 4 ) by using the elementary transformations. To be explicit, we list the 
relevant relations 


Y u (B. 2 )=V 23 Y u (B 1 ) 


Y u (B 3 ) = V 12 V 23 Y v (B 1 ) 


n(B i+3 ) = n(B i )Q 


( 10 ) 


where the index i = 1,2,3 in the last equality is implied. The same transformations apply to the patterns 
in Eq. Q. As we will show later in this section, the above observations will be useful to analyze the texture 
zeros in the effective neutrino mass matrix M u . Note that the elementary transformations are implemented 
to examine the location of texture zeros, so the nonzero elements in both Y v and the corresponding M v 
are not necessarily identical for each pattern. 

Below the seesaw scale, one can integrate out heavy Majorana neutrinos and obtain the unique Wein¬ 
berg operator 0 3 = (/c/2) (£ L H) ■ (H T ££) of dimension five [49j with k = —Y u M^ x YJ. After the spon¬ 
taneous gauge symmetry breaking, neutrinos acquire tiny Majorana masses from the Weinberg operator 
and their mass matrix is M u = nv 2 , which is just the seesaw formula in the language of effective theories. 
Now it is clear that Y v is given at a superhigh-energy scale /r = M i , but neutrino oscillation parameters 
are measured at low energies. In order to study whether the flavor structure of Y u in Case A, B and C 
is viable, we have to examine the RG evolution of k from the seesaw scale A gg to the electroweak scale 
A ew , and compare the predictions from k(A ew ) with neutrino oscillation data. 

Given Y v in Eqs. (0) ((9|) , we are ready to check if k inherits some texture zeros from Y v . Since all the 
patterns in each class are related by V ij and Q matrices, it is sufficient to consider the first pattern and 
perform the corresponding elementary transformations to derive the results for the others. More explicitly, 
we have k{M x ) at the seesaw scale 


( 11 ) 


where one can observe that the patterns C i (for z = 1,2,--- , 6 ) do not lead to any texture zeros in k. For 
Case A in Eq. (0), it is easy to derive k a = 'P l jK A V 1] for i = 2,3, so k in this case has a nonzero 2x2 


: 

/0 0 o\ 

0 x x 

’ Kb i : 

fx 0 x\ 

0 x x 

> : 

(X X x\ 

XXX 


X 

X 

0 


^X X Xj 


\X X X) 
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block submatrix. For Case B in Eq. ([8jh with the help of Eq. m, we arrive at the following identities 


k B 2 — ^*23^^23 i k B 3 — ^ , 12^ , 23 k B 1 ^ , 23^ , 12 ) K B i+3 “ K B ; ) (12) 

where the last identity indicates that one texture zero is located in the same position in k for B- +3 and 
Bj- for i = 1, 2,3. 


2.3 Renormalization-Group Running Effects 


As we have mentioned, neutrino masses at the sub-eV level indicate that the seesaw scale is extremely 
high A ss 10 14 GeV, if the Dirac neutrino Yukawa couplings are of order 0(1). In the full theory above 
the seesaw scale, two heavy Majorana neutrinos are added into the SM particle content, and they interact 
with the SM particles only through the Yukawa interaction, which is governed by the coupling matrix Y u . 
After taking into account radiative corrections and renormalizing the model in the scheme of dimensional 
regularization and modified minimal subtraction, we are left with coupling and mass parameters that 
depend on the renormalization scale y. The evolution of model parameters with respect to /i is described 
by their RG equations. For y < A gs , the decoupling of heavy Majorana neutrinos is treated by explicitly 
integrating them out, and the low-energy effective theory turns out to be just the SM plus a dimension- 
five operator, which is responsible for neutrino masses. At the one-loop level, the RG running effects of 
neutrino masses and flavor mixing parameters can be studied by solving the RG equation of k [501 EH [52] 



with t = ln(///A EW ). In the SM, the relevant coefficients in Eq. (fl3l) are C K = —3/2 and a K ~ —3 g\ + 
6 y/ + A, where g 2 stands for the SU(2) L gauge coupling, y t the top-quark Yukawa coupling, and A the Higgs 
self-coupling constant. If the dimension-five Weinberg operator is derived in the minimal supersymmetric 
standard model (MSSM), we have M v = K(usin/3) 2 with tan/3 being the ratio of vacuum expectation 
values of two MSSM Higgs doublets. In this framework, the RG equation of k is still given by Eq. m 
but with C K = 1 and a K ~ —6gf/5 — 6 g 2 + 6y/. Note that only the top-quark Yukawa coupling is retained 
in a K , as the Yukawa couplings of other fermions are much smaller and have safely been neglected. The 
RG evolution of neutrino masses and lepton flavor mixing parameters has been extensively studied in the 
literature Ell E51 ESSES EH- See, e.g., Ref. [58], for a recent review on this topic. 

Working in the basis where the charged-lepton Yukawa coupling matrix Y l = diag{y e , y^, y T } is diagonal, 
we can solve Eq. (fT3l) and obtain 


k (Aew) — h 

(i e 0 o\ 
0 /„ 0 


(l e 0 o\ 
0/,0 


G o ij 


" 0 ij 


(14) 


where the evolution functions read 

I 0 = exp 

I a = ex P 


1 

167T 2 


167T 2 


L 

l 


ln(Afj/A 

In(Mj/A 


Ew) 


Ew) 


a K (t) df , 

2/a(t) d t , 


(15) 

(16) 


for a = e, y, r. From Eq. it is now evident how the low-energy observables residing in M u = k(A ew )v 2 
are related to the model parameters in k(M 1 ) at a high-energy scale. In the following, we show that it is 
already possible to exclude most patterns in Eqs. (J7])- ([9j) based on the solution in Eq. (ITT]) . 












1. An important observation from Eq. (1141) is that texture zeros in k are rather stable against the RG 

running. On the other hand, Eq. dill) tells us that k(M 1 ) for the patterns A - possesses five vanishing 
elements, appearing in the i-th row and i-th column. Therefore, k(A ew ) in Case A inherits the same 
structure of leading to just one nontrivial mixing angle, which has already been excluded by 

current neutrino oscillation data. Thus, all three patterns in Eq. 0 are ruled out. 

2. Then we turn to the patterns B x 2 3 , and the same conclusions should also be applicable to B 4 5 6 , 
since the texture zero in M u is located in the same position. For this class, there is only one texture 
zero in k(A ew ) or M v = k(A ew )u 2 in the off-diagonal position, namely, 

( M u)a/3 = m i U *i U Pi = 0 > ( 17 ) 

i 

for (a, (3) = (e,/i), (e,r) and (n,r). When the RG running effects are considered, Eq. (fTTl) indicates 
that the texture zero remains in the effective neutrino mass matrix M u . The constraints on neutrino 
masses and mixing matrix elements in Eq. ED can be expressed as 


U a2 U B2 m 2 + U a3 Up 3 m 3 = 0 for NO 

U al U f31 m l + U a2 U 32 m 2 = 0 for 10 


(18) 


which have been investigated in Ref. m , where the latest neutrino oscillation data are implemented 
but the RG running effects are entirely ignored. In the NO case, it has been found that all the 
patterns in Eq. (j 8 j) are ruled out mainly due to the observed 9 13 [59J [60, .EE E2]- In the 10 case, 
(M u )^ r = 0 is shown to be strongly disfavored, so the patterns B 3 and B 6 are excluded. Hence, 
according to Ref. m, only B 3 2 and B 15 in the 10 case are compatible with the latest neutrino 
oscillation data. 


3. Since the patterns in Eq. @ do not imply any zero elements in k(M 1 ), the analysis of Case C 
in Ref. JT3] seems to be not applicable. Thus it is expected the predictions at a superhigh-energy 
scale will be significantly changed at the low-energy scale. However, as we demonstrate below, 
a characteristic relationship among the elements in k is maintained at the low-energy scale and 
validates the conclusions in Ref. [T3j. Let us take the pattern C 3 for example, and specify its matrix 
elements: 


Cr 



k(M x ) 


1 

mJ" 


/o 0 0 \ 
0 0 0 
\0 0 a' 2 ) 



7 a 2 ab ab'\ 
ab b 2 W 
y,ab ' bb' b ,2 J 


(19) 


where the corresponding k(M 4 ) has been given as well. Combining Eq. (1141) and Eq. (1191) . one can 
verify that the relation 

(M„) ee = = (M„) er 

(^),e (M^ r 1 j 

holds both for [i = A EW and for fj, = M 1 . Therefore, it is adequate to inspect if the relationship 
in Eq. (1201) is satisfied by current neutrino oscillation data. More explicitly, the first identity in 
Eq. (1201) gives rise to U e3 U /l2 = U e2 U^ 3 for NO, and U e2 U. j3 = L r el t / /i2 for IO, while the second 
identity is fulfilled automatically. The constraints for the other patterns can be found in a similar 
way. Those relations among the PMNS matrix elements have also been derived in Ref. [13], although 
in a different manner, and used to exclude all the patterns in Eq. Q in both NO and IO cases. 
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In summary, we have proved that texture zeros or proportionality relations in k(M 1 ) are not spoiled 
by the RG running effects, so they also exist in k(A ew ) at the low-energy scale. Consequently, neutrino 
oscillation data can be directly implemented to rule out most patterns of Y u with two texture zeros. It 
turns out that only B 4 2 and B 4 5 in Eq. ([8]) in the case of 10 are consistent with experimental data, which 
generalizes the conclusions reached in Ref. m to the situation including radiative corrections. 


2.4 Viable Patterns 

Now we are left with just four viable patterns, namely B 4 2 and B 45 in Eq. ([8]), and only the 10 case is 
allowed. The latter indicates a sizable value of m^, around 50 meV, and thus is quite encouraging for 
future experiments to search for neutrinoless double-beta decays. Although the RG running effects are 
unable to revive any patterns in the NO case, they do have significant impact on the allowed regions of 
model parameters, particularly in the MSSM with a large tan/3. Hence, in this subsection, we examine 
four viable patterns in more detail, and explore the favored parameter space. 

As we have shown in the previous subsections, the effective neutrino mass matrix M u at the low-energy 
scale in this case contains one texture zero, which sets two constraining relations on neutrino masses and 
mixing angles. Since neutrino mass spectrum is completely fixed by the observed neutrino mass-squared 
differences, one can determine two CP-violating phases in terms of neutrino masses and three mixing 
angles. According to Eq. m , the two patterns in each pair of {B 1; B 4 } and {B 2 ,B 5 } are related by 
an exchange between two columns, so the location of texture zero in M u is identical, indicating the same 
low-energy predictions. However, the model parameters in the full theory at the seesaw scale are different, 
as we shall show later. Using the second identity in Eq. (fT8l) for the case of (a,/3) = (e,//), we obtain 

m l c 12( c 23 s 12 + c 12 s 23 s 13 el<S ) _ m 2 s 12( c 12 c 23 — s 12 s 23 s 13 el<5 )e 2l ' T = 0 > (21) 


whose real and imaginary parts allow us to determine 5 and a via 


cos 5 = 

cos 2cr = 


444(i-4+ 44(4-44) 

2' s 12 C 12' S 23 C 23 S 13(4 + 44) 

444(i + c 2 )- 44(4 + c 2 4) 

2 C44(4 +44) 


( 22 ) 

(23) 


up to a sign ambiguity. Since 1 — £ 2 ~ \/2s 2 3 ~ 0.03 holds as an excellent approximation, one can expand 
the right-hand sides of Eqs. (12211 and (12311 in terms oil — C, 2 and s 2 3 , and ignore the higher-order terms of 
0(s\ 3 ). After a straightforward calculation, we arrive at 


cos 5 

cos 2(7 


sin2 4 ^2 

4 tan @ 23 sin# 13 

tan 2 d 23 sin 2 0 l3 
2 sin 2 0 12 cos 2 d 12 


tan 0 23 
tan 20 12 


sin # 13 , 


(24) 


implying that 5 ~ 90° or 270° and <7 ~ 0°. The deviation of <5 from the maximum 90° or 270°, and that 
of c from zero, are on the order of 0 l3 in the leading-order approximation. For the pattern B 2 , one needs 
to consider Eq. (fT8l) with (a,/3) = (e, r). It is easy to verify that Eqs. (1221) (1+41) become applicable to this 
case after replacing 5 with <5 + 7r, as well as d 23 with ir/2 — 0 23 , namely, flipping the octant of 0 23 . This 
observation indicates that the determination of the octant of d 23 and the measurement of CP-violating 
phases 5 in future neutrino oscillation experiments can be used to distinguish between the patterns B 4 
and B 2 for the Dirac neutrino Yukawa coupling matrix. 
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There are five real parameters in M D , since two matrix elements are zero and three arbitrary phases 
can be absorbed by redefining the charged-lepton fields. Moreover, the heavy Majorana neutrino masses 
M 1 and M 2 are free parameters. It is convenient to introduce the Casas-Ibarra parametrization [63] 


M d = U\! M u O\j M r = U 



0 ^ 
0 


/ cos 2 — sin z'' 
sin 2 cos z 


V 0 0 / 



(25) 


where U is the PMNS matrix given in Eq. ©> and O is a 3 x 2 orthogonal matrix with z being a complex 
parameter, satisfying O t O = OO 3 = 1. Note that we have concentrated on the 10 case, which is the only 
allowed possibility in the FGY model. All the mixing angles, CP-violating phases, and neutrino masses 
in Eq. (|25|) should take values at the seesaw scale, which are in general distinct from those extracted from 
neutrino oscillation experiments at the low-energy scale (e.g., at the Fermi scale M z = 91.2 GeV). Because 
of the texture zeros in M D , the CP-violating phases 5 and a can be determined in terms of neutrino masses 
and mixing angles as in Eqs. (1221) and (l23lh but now with their values at the seesaw scale. In addition, the 
complex parameter z can be determined by 


tan z = 



tan 0 12 


(26) 


for B 1 and B 2 . Since B 4 and B 5 are related to B 4 and B 2 by exchanging two columns, respectively, the 
parameter z in the former two cases can be calculated first from Eq. (I26jh and then followed by a shift 
of z — y z T 7t/2. Now it is evident that the complex parameter z is actually determined by the neutrino 
mass ratio £ = m 1 /m 2 , the mixing angle d l2 and the Majorana CP-violating phase a. However, the RG 
running effects on these parameters, in particular 6 12 and er, could be significant. 

Taking Pattern B x for example, we proceed to explore the possible parameter space at the low-energy 
scale by using the global-fit results in Table 1, and that at the high-energy scale by numerically solving the 
complete set of one-loop RG equations. In view of minimality of the FGY ansatz, we shall consider the 
minimal SM. In the SM, the largest charged-lepton Yukawa coupling y T is as small as 10 -2 . According to 
Eq. (1161) . the evolution function running from the electroweak scale to the seesaw scale A gg = 10 13 GeV is 
approximately given by I T ~ exp(—25 x ICG 6 ) ~ 1. Therefore, we have I e ~ I fL ~ I T ~ 1, and the form of 
n remains unchanged during the RG running, resulting in negligible modifications on the mixing angles, 
CP-violating phases, and the ratio of neutrino masses. This means that the predictions of FGY ansatz 
are essentially valid at high-energy scales in the minimal SM. 

In the MSSM, the running effects are expected to be significant, since y T can be enhanced by large 
values of tan/3. We first input the neutrino mixing angles and two neutrino mass-squared differences 
within their 3<r ranges at M z . Two stages of RG running are then performed, namely, one from M z to the 
sparticle mass scale Mgugy with the SM RG equations, and the other one from Mgugy to A gg = 10 13 GeV 
by adopting the MSSM RG equations. Taking MgusY = 1 TeV, we have calculated the running effects on 
neutrino mixing parameters, and the numerical results are presented in Fig. [I] We have also tried to vary 
this intermediate sparticle mass scale MgusY from 1 TeV to 10 TeV, however, only minor changes (< 5%) 
are found on the mixing parameters. 

In Fig. [T] the allowed regions of three neutrino mixing angles {0 12 , $ 13 ) @ 23 }; two leptonic CP-violating 
phases {<5,< 7 } and two nonzero neutrino masses {m 1 ,m 2 } are shown in the MSSM with tan/3 = 30 and 
tan/3 = 50. The allowed parameter space at the low-energy scale is denoted by black points, and one can 
observe that <5 and a are restricted to a small area around b = 90° and a = 10°. This observation can 
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Figure 1: Illustration for the RG running effects on neutrino mixing angles {#121 $131 $ 23 }) leptonic CP- 
violating phases {<5, <r} and neutrino masses {m 1 ,m 2 } for Pattern B x in the MSSM, where the black 
points denote the parameters at M z = 91.2 GeV, while the dark- and light-gray points represent the 
parameters at the seesaw scale A gg = 10 13 GeV for tan /3 = 30 and tan/I = 50, respectively. Note that 5 
and <7 also have another branch of solutions with their signs inverted simultaneously, and the mass scale 
of sparticles is taken to be M gugY = 1 TeV. 


be easily understood with the help of Eq. (1241) . which indicates that the deviations of (d, a) from (90°, 0°) 
are measured by the neutrino mass-squared difference Am^ = (1 — ( 2 )m 2 and the small but nonzero 
mixing angle 0 13 . At the high-energy seesaw scale A gg = 10 13 GeV, the parameter space in the MSSM 
with tan/3 = 30 and tan/3 = 50 has been represented by dark- and light-gray points, respectively. One 
can see that the RG running effects on # 13 and 0 2 3 are insignificant, whereas the running effects on 0 12 , 5 
and a are indeed remarkable. Therefore, it is necessary to include the running effects on those parameters 
when we consider the generation of baryon number asymmetry in our Universe, which takes place at a 
superhigh-energy scale. 

From Eq. (1261) . we can figure out the real and imaginary parts of z in terms of neutrino mixing 
parameters. More explicitly, we have 


R“ K “2 

Imz ~-In 

4 


arctan 


sin a + cot 6 


12 


— arctan 


cos a 

1 — 2 sin a cot 0 12 + cot 2 9 


sin a — cot 0 l2 


cos a 


12 


(27) 


1 + 2 sin a cot 0 12 + cot 2 0 12 j 

where £ ~ 1 is assumed. For a small tan /3, the RG running effects are negligible, so the mixing parameters 
can be identified with those extracted from oscillation experiments. In this case, one can expand Eq. (127|) 
in terms of the Majorana CP-violating phase a, which is constrained to be small. At the leading order, 
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Figure 2: The absolute value \z\ and the phase argz of the complex parameter z are given in units of 
degrees in the left and right panels, respectively. For a given value of tan/3, \z\ and arg z are calculated 
by varying the low-energy parameters in their 3tr ranges and the high-energy scales from 10 8 GeV to 
10 13 GeV. Note that 2 is almost real in all cases as indicated in the right panel. 


we get \z\ ~ 7r/2 — 0 12 and aigz ~ <rsin2(9 12 /(20 12 — n). In the other extreme case, where the RG running 
is significant for a large tan/3, we can expand Eq. (1271) in terms of 0 l2 and obtain \z\ ~ n/2 — cos a tan d 12 
and arg z = 2sincrtan0 12 /7r. In both cases, arg 2 is found to be close to the real axis, i.e., around 5°. In 
general, both o and 0 12 are not small angles, and the above approximations are invalid. 

However, one can compute the complex parameter z by inputting the low-energy values of neutrino 
mixing parameters and solving the RG equations. The numerical results of \z\ and arg z are given in 
Fig. El where both small and large values of tan/3 are considered. Furthermore, the 3a ranges of mixing 
parameters and a seesaw scale within [10 8 ,10 13 ] GeV are taken into account. One can see from the right 
panel of Fig. [2] that a small phase of z is obtained in all cases, implying the suppression of CP violation 
at the high-energy scale. The latter observation becomes clearer when we calculate the CP asymmetries 
in the decays of heavy Major ana neutrinos. 


3 Baryon Number Asymmetry 

One salient feature of the canonical seesaw model is to simultaneously explain tiny neutrino masses and 
the observed baryon number asymmetry in our Universe, which is usually measured by the baryon to 
photon density ratio [S3] 

r/g = —^ = (6.065 ± 0.090) x 10" 10 , (28) 

n 7 

where n B and n 7 stand for today’s baryon and photon number density, respectively. In the very early 
Universe, when the reheating temperature after inflation is so high that heavy Majorana neutrinos N i can 
be produced in thermal equilibrium. As the Universe cools down, the CP-violating decays of N i will go 
out of thermal equilibrium if the decay rate becomes smaller than the expansion rate of the Universe. The 
CP asymmetries in the decays of N i into leptons of different flavors are defined as [271 EH EH 

= r(y t ^z Q g)-r(A^7 Q i?) 

F(A7 -> l a H) + T(A7 l a H) ’ 

where T(N t —>• l a H ) and r(A^ —>• l a H) for a = e,fi,T denote the decay rates of N i into leptons l a and 
anti-leptons l a , respectively. It is the interference between the tree-level and one-loop decay amplitudes 
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that gives rise to CP asymmetries, which receive both contributions from the one-loop self-energy and 
vertex corrections. More explicitly, we obtain 


&ia 


1 

87 T(YjY v )u 


im^TOO^u 

k^i 



(30) 


where x ki = M k /Mf and the loop functions are defined as follows 


f( X ki) 

9(x ki ) 


'%ki 


1 — X 


ki 


(! - x ki ) 2 + r2 


+ 1 


ki 


1 — X 


ki 


(! " x ki ) 2 + C 2 


ki 


(1 + x ki )ln 


1 + %ki 
%ki 


(31) 


If the mass spectrum of heavy Majorana neutrinos is strongly hierarchical, r ki can be neglected in the 
denominators in Eq. m- However, it serves as an important regulator to avoid any singularity in the 
limit of mass degeneracy M? = Mf or equivalently x, - = 1. In the resonant regime, the true form of r ki 
is still controversial at present [65], and three distinct expressions have been derived: (i) r ki = x ki T k/ M k 
by a quantum field-theoretic approach j66[ 1ST] : (ii) r ki = Y i /M i — x ki T k /M k by a modified version [681 69j 
of the approach introduced in Ref. [66]; (iii) r ki = T i /M i + x ki T k /M k by an effective Kadanoff-Baym 
approach with a specific quasi-particle ansatz mm As we numerically demonstrate in the FGY model, 
three different expressions of r ki lead to the same result if a successful leptogenesis is realized. 

The produced lepton-number asymmetries in the N i decays will partly be washed out by the inverse 
decays and lepton-number-violating scattering, if these processes proceed efficiently. In order to describe 
the washout effects, we introduce the decay parameters = T i /H(M i ), where Tj = {YyY v ) ii M i /S'K is the 
total decay width of N i and H(M i ) is the Hubble parameter at temperature T = M i . In the radiation- 
dominated epoch, the Hubble parameter is given as a function of temperature H(T ) = 1.66-^/$* (T)T 2 /M pl , 
where M pl = 1.2 x 10 19 GeV is the Planck mass and g*(T) is the number of relativistic degrees of freedom 
at T. The lepton number asymmetries will be converted into the baryon number asymmetry through the 
(B + L)-violating and (B — L)-conserving sphaleron processes [25] [26], which are in thermal equilibrium 
between T = 200 GeV and 10 12 GeV. The final baryon number asymmetry is then given by m 


r] B » -0.96 x 10 2 ( 32 ) 

i OL 

where the efficiency factors n ia can be determined by solving the Boltzmann equations of heavy Majorana 
neutrino and lepton number densities. Roughly speaking, they are governed by the flavor-dependent decay 
parameters K ia = P ia K i , where P ia = \(Y u ) a i \ 2 stands for the projection probability of the final 
lepton state in N { decays onto a specific lepton-ffavor state. 

So far, we have focused on leptogenesis in the SM. In the MSSM, the CP asymmetries in the decays 
of both A] and its superpartner are twice larger, since the number of particles running in the loops are 
doubled. However, in the strong washout regime, the inverse decay rates are also doubly efficient, reducing 
the lepton asymmetries by a factor of two. In addition, the particle content is twice much in the MSSM, 
so we have the number of relativistic degrees of freedom g* = 228.75 in the MSSM, while g* = 106.75 
in the SM. Altogether, the baryon number asymmetry in either strong or weak washout regime in the 
supersymmetric case is not much changed with respect to the non-supersymmetric case [29] . 

In the vanilla scenario of leptogenesis, the mass spectrum of heavy Majorana neutrinos is taken to be 
hierarchical, and only the lightest Majorana neutrino A r 1 and the one-flavor approximation are considered. 
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This is actually done for the FGY model in the previous papers [201 [2TL 22i 113j . where a narrow mass 
range of the lightest heavy Majorana neutrino M 1 ~5 x 10 13 GeV has been found in the 10 case. In the 
following, we calculate the baryon asymmetry via a flavor-dependent leptogenesis by taking into account 
the lepton flavor effects and non-hierarchical mass spectra of heavy Majorana neutrinos. 


3.1 Lepton Flavor Effects 


The interaction rates associated with charged-lepton Yukawa couplings become larger than the expansion 
rate of the Universe at different temperatures, and thus affect the washout effects on lepton number 
asymmetries [7211731 El [751 [76]. For > 10 12 GeV, the leptogenesis mechanism works at the temperature 
T ~ where all the charged-lepton Yukawa interactions are negligible compared to the expansion 
rate. Therefore, the lepton state produced in the decays also participates in the inverse decays and 
lepton-number-violating scattering. In this case, it is valid to treat leptons as a single flavor in both 
generation and washout of lepton number asymmetries. The relevant quantities are just the total CP 
asymmetry = Yl a £ ia and the efficiency factor k which is determined by the decay parameter K i . 
For 10 12 GeV > M i > 10 9 GeV, the r charged-lepton Yukawa interaction is in thermal equilibrium and 
able to single out the r lepton flavor in the thermal bath. Therefore, one has to deal with two lepton 
flavors, namely the r flavor and a combination of e and /i flavors. The relevant parameters are the CP 
asymmetries s iT and s i2 = £ ie +£ ifX , and the efficiency factors K ir and n i2 , which are calculable by using K ir 
and K i2 = K ie + K ifl . For M x < 10 9 GeV, both r and /i charged-lepton Yukawa interactions are efficient 
enough to recognize r and // flavors in the system, implying that a three-flavor treatment is necessary. 

First, we compute the CP asymmetries in the FGY model. Since the Dirac neutrino Yukawa coupling 
matrix is given in Eq. (1251) . it is straightforward to figure out £ ia in Eq. (1301) . In the hierarchical limit 
of Afj -C M 2 , we need to just focus on e lQ and assume that the lepton asymmetries generated from the 
decays of N 2 have been washed out by the ^-related lepton-number-violating processes. For Pattern B 1 
with (Y v ) el = {Y v ) ij2 = 0, we obtain e le = £ lfl = 0, and 


3 M l Am^Iinfc 2 ] 
16-7T v 2 m^cj 2 + m 2 \s z 


(33) 


where the second equality has also been found in Ref. m ■ The CP asymmetry is suppressed by the tiny 
neutrino mass-squared difference A m 21 ~ 7.5 x 10 -5 eV 2 . Furthermore, as we have shown in the previous 
section, the complex parameter z is very close to the real axis, implying that |Im[c 2 ]| ~ |^| sin(2|z|) arg(z) 
should also be small. The numerical values of |Im[c 2 ]| have been presented in Fig. [3] for a wide range of 
model parameters, where one can observe that |Im[c 2 ]| is actually small and varies between 0.03 and 0.09. 
In the present work, we shall concentrate on Pattern B 1; but one can calculate the CP asymmetries for 
the other three viable patterns in a similar way. The important results for all four viable patterns have 
been summarized in Table 2. 

Second, instead of solving the complete set of Boltzmann equations, we apply the analytical formulas 
obtained in Ref. m to estimate the efficiency factors. If the initial thermal abundance of heavy Majorana 
neutrinos is assumed, the efficiency factor is approximately given by m 




1 — exp — 




(34) 


where z B (K ia ) = 2 + 4AT? Q 13 exp(— 2.5/K ia ). Hence the efficiency factors are completely fixed by the decay 
parameters K ia , which are in turn determined by the flavor structure of Y u . For Pattern B x , we get the 
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Figure 3: Numerical results of Im[c 2 ] are calculated by solving the RG equations with the low-energy 
parameters in their 3 <t ranges as inputs. The high-energy scale has been chosen to be 10 8 GeV (red stars), 
10 10 GeV (green dots), and 10 13 GeV (blue triangles) for a given tan/3 in the MSSM. 


total decay parameter 


K i = 


M 2 v 2 (m 1 \c z \ 2 + m 2 \s z \ 2 ) 


50 , 


(35) 


where Imz <C 1 and m 2 ~ m 1 
determined by 

P\r 


0.05 eV have been used in the last step. The projection probability is 


lOQ-rll 2 = \ U rlV^- C z + Ur2V^ S z\ 2 = . n 
\( Y u)^i \ 2 \U^y/rnlc z + U^ 2 y/mis z \ 2 ' 


(36) 


where the identity tan z = —U ely /mi/(U e2y /m 2 ) has been implemented to significantly simplify the result. 
Given d 23 ~ 45°, we arrive at P 1t ~ 0.5 and K lr ~ 25. For comparison, we can also figure out P lr for 
Pattern B 4 with (Y v ) e2 = {Y v )^i = 0. With the constraint tanz = —U fjls /mi/(U ll2 y/m 2 ), we have 


Pir = l( y Jril 2 = \ U riVPnc z + U T 2 y /mjs z \ 2 _ tan 2 6> 13 
1 - P\t l(^)el| 2 l^el\MTc 2 + U e 2 y/m^S z \ 2 ~ COs6» 23 

and thus P 1t ~ 0.05 and K lr = 2.5, which are one order of magnitude smaller than the result in the 
previous case. Since P lr in Eq. (1361) or Eq. (1371) depends mainly on 0 23 and 0 13 , its value should be quite 
stable against the RG running. 

With both the CP asymmetries and decay parameters, we are ready to find out the efficiency factors, 
and then baryon number asymmetry. The numerical results are summarized as follows: 


• If M 1 > 10 12 GeV, we can treat leptons as a single flavor, and the relevant quantities are the CP 
asymmetry ~ —2 x 10 -6 (M 4 /10 13 GeV), which is identical to e lT as shown in Eq. (1331) . and the 
efficiency factor k 3 ~ 5 x 10^ 3 by inserting K x = 50 into the analytical formula in Eq. (1341) . Putting 
all together, we obtain the baryon number asymmetry 

Vb « -0.96 x 10" 2 £ 1 k 1 = 1.0 x 10” 10 ( 1Ql fg eV ) , (38) 

which is in agreement with the result in Ref. m- Therefore, heavy Majorana neutrinos should be as 
heavy as 6 x 10 13 GeV to generate the correct baryon number asymmetry. Note that |Im[c 2 ]| = 0.05 
has been assumed in the above calculation, but it is evident from Fig. [3] that the RG running effects 
on mixing parameters can enhance or reduce this value by a factor of two, depending on tan /3. 
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• If Mj < 10 12 GeV, the CP asymmetry is given by the same formula e lr « —2 x 10 -7 (M 1 /10 12 GeV), 
which will be at least one order of magnitude smaller compared to the previous case. Since the flavor 
structure of Y v under consideration indicates s le = £ lfl = 0, there is no contribution from other 
lepton flavors to the lepton number asymmetries. The washout of lepton number asymmetries is 
now determined by K 1t = P lr K 1 = 25, leading to an efficiency factor k 1t ~ 0.01. Although there is 
an enhancement by a factor of two, the mass of the heavy Majorana neutrino is too small to provide 
a large enough CP asymmetry. If we turn to the case of Pattern B 4 , the CP asymmetry remains 
the same and the efficiency factor is k 1t ~ 0.2, so we have the final baryon number asymmetry 

Vb ~ -0-96 x l(T 2 e lr K lT = 3.8 x 1(T 10 , (39) 

which is on the right order of magnitude even for M 1 = 10 12 GeV. However, it is worthwhile to 
point out that M 1 = 10 12 GeV is on the edge of two-flavor approximation, when the coherence of 
lepton state in N 1 decays may be destroyed by the r Yukawa interaction. In this case, the classical 
Boltzmann equations are not accurate enough to give the correct answer, and the fully quantum 
Boltzmann equations should be applied isa izn Eg. Hence the flavor effects may open a possibility 
to realize a successful leptogenesis even for a smaller M 1 . 

For even smaller masses M 1 <C 10 12 GeV, the CP asymmetries are significantly suppressed. It is impossible 
to explain the observed baryon number asymmetry in the FGY model, although the flavor effects tend to 
protect lepton number asymmetry from washout. 


3.2 Beyond Hierarchical Limit 


The high mass scale of heavy Majorana neutrinos causes the so-called naturalness or fine-tuning problem 
for the light Higgs boson mass ES EH EMEU, and the gravitino overproduction problem if the model 
is supersymmetrized [35]. In Ref. ITT , a detailed analysis of the naturalness problem in the type-I seesaw 
model yields an upper bound on the heavy Majorana neutrino masses, namely, M 1 < 4 x 10 7 GeV and 
M 2 < 7 x 10' GeV. These upper bounds have been derived by requiring that the radiative corrections 
induced by heavy Majorana neutrinos to the Higgs boson mass should be around the TeV scale. Obviously, 
this bound is in contradiction with the requirement of M 1 ~ 10 13 GeV for explaining the baryon number 
asymmetry in the FGY model. Therefore, it is well motivated to go beyond the hierarchical limit and 
consider both mild mass hierarchy and a nearly-degenerate mass spectrum. 

In the mild hierarchy case, we take M 2 to be a few times M 1 . For the later convenience of quantifying 
the level of mass degeneracy, we introduce a dimensionless parameter 


A = 


M 2 - M 1 

m 2 


(40) 


which is zero in the limit of exact mass degeneracy M 1 = M 2 and approaches one for M 2 3> M 1 , which is 
the case discussed in the previous subsection. 

Because of a mild hierarchy between and M 2 , both N 1 and N 2 participate in the production and 
washout processes of lepton number asymmetries. The evolution of these asymmetries therefore involves 
solving the Boltzmann equations with both N 1 and N 2 , and the previously used analytic formula for 
estimating the efficient factor is no longer applicable. To obtain a rough estimation of the baryon number 
asymmetry in this mild hierarchy case, we next consider a simplified set of Boltzmann equations, where 
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only the inverse-decay processes are included in the washout term. First, the evolution equations of N-± 
and N 2 number densities are m 


dn 


N ; 


dz 


= ~ D i{n Ni ~n$) 


(41) 


where z = M x /T, and n N . is the number density for N t normalized by its density in ultra-relativistic 
thermal equilibrium (i.e., T S> MB. Here n e 3 = zfK,2{zi)/2 with z% = MJT = zMi/Mi is the density in 
thermal equilibrium, and /C 2 (z) is the modified Bessel function of the second kind. The decay factor D{ is 
defined to be 


H {2 


(42) 


where K i has the same form as the previously defined total washout factor, and (l/yj = IC 1 (z i )/IC 2 (z i ) 
is the thermally averaged dilation factor. Second, we also have the evolution equations for the lepton 
asymmetries, namely, 



£' 


x Di(n N - n eq 


N. 


n A Q E 
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la " 1 


ID 


(43) 


where n Aa is the B — L asymmetry density for the flavor a, which has also been normalized by the density 
of Nj in the ultra-relativistic thermal equilibrium, and the total B — L asymmetry density n B _ L is then 
given by n B _ L = n Aa . In addition, Pi a is the projection probability defined previously, and the 
inverse-decay washout term VF/ d is as follows 


W} D = . 

4 Mi 


(44) 


Given the above set of Boltzmann equations, we then solve them numerically. The initial conditions 
are obtained by setting the thermal abundance of n N , and vanishing B — L asymmetries. In Fig. U we 
present the allowed parameter space for M 1 and A in the case of Pattern B 1 . The black solid curve 
represents a contour of r] B = 6.065 x 10~ 10 , for which the observational uncertainty is so small that it will 
be hidden by the line width in the figure. The mass regions, which are represented by the shading areas, 
are characterized by the charged-lepton flavor effects. 

In the highly degenerate case, we calculate rj B in two ways: solving the simplified set of Boltzmann 
equations introduced Eqs. (Hll) and (l43l) . and applying the approximate analytical formulas. In [78] , it 
was argued that in the degenerate limit, the N 1 and N 2 washout contributions add up, resulting in 


7 7b — —0.96 x 10 2 + £ 2 a) K (Kia + K 2a ), (45) 

OL 

where the efficiency factor k is still calculated via Eq. (1341) . The summation over a depends on the region 
of the lepton flavor effects. We focus on Pattern B 1 with (Y^) e i = ( Y u ) 2 = 0, and the other cases can 
be analyzed in a similar way. 

In Fig. [J1 we show the allowed parameter space for M 1 and A for a variety of masses, in the single- 
flavor, two-flavor and three-flavor regions. We have demonstrated that the two approaches with simplified 
Boltzmann equations and approximate formulas lead to the same result. In the mild hierarchy case, we 
observe from Fig. 0] that M l still sits around 5 x 10 13 GeV. This can be easily understood, as we know 
that e 2 is at most as large as e 1 . To see this point clearly, we first calculate e l /e 2 by using Eq. (|29l) . and 
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M : (GeV) 


Figure 4: Illustration for the dependence of baryon number asymmetry on the lightest heavy Majorana 
neutrino mass M 1 and the mass degeneracy parameter A. The black and solid curve corresponds to the 
allowed regions of model parameters, for which the observed baryon number asymmetry r/ B ~ 6.065 x IIP 10 
can be naturally explained. The dashed lines indicate a few typical values of the mass ratio M 2 /M 1 . 

find it divergent when e 2 = 0, corresponding to M 2 /M 1 ~ 2.36. When M 2 /M 1 < 2.36, e 1 and e 2 have 
the same sign, while the opposite situation happens when M 2 /M 1 > 2.36. In addition, |e 1 | > |e 2 | holds 
for all ratios of M 2 /M 1 . Therefore, including the contributions from N 2 cannot significantly enhance the 
amount of CP asymmetry, and one then still needs to raise the mass scale of M 1 so as to reach the required 
value of ? 7 b . In the nearly-degenerate case, we see that a mass degeneracy at the level of A = 10~' is 
required to meet the naturalness bound M 1 < 4 x 10' GeV and account for the baryon number asymmetry 
via resonant leptogenesis [66l ’67, [79]. In our calculations, the formulas of CP asymmetries with different 
regulators r ik lead to the same result in the FGY model. Although it seems unnatural to require such a 
high mass degeneracy, it can actually be achieved by implementing a flavor symmetry and its soft breaking 
at a superhigh-energy scale m , or by the RG running effects [80, IHT] . As one can see, there is a kink 
around M 1 = 10 12 GeV. The reason is simply that we use different Boltzmann equations for the two 
cases of below and above 10 12 GeV. The kink should disappear if the fully quantum Boltzmann equations 
with coherent flavor effects are used [65]. The curve is continuous around M 1 = 10 9 GeV, since the flavor 
structure enforces only one nonzero CP asymmetry e lr . 

4 Summary 

In light of the latest neutrino oscillation data, we have performed a further study of the FGY model, in 
which only two right-handed heavy Majorana neutrinos are introduced and two texture zeros appear in 
the Dirac neutrino Yukawa coupling matrix, by taking into account the RG running of neutrino mixing 
parameters and flavor effects in leptogenesis. Such an investigation is well motivated in two aspects. 

First, the FGY model is very interesting and predictive, and can be readily confronted with the latest 
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neutrino oscillation data. Since the lightest neutrino is massless, the neutrino mass spectrum is fixed by the 
neutrino mass-squared differences, which are precisely measured in neutrino oscillation experiments. There 
are one Dirac and one Majorana CP-violating phases, which are actually determined by neutrino mixing 
angles and masses. The neutrino mass ordering is inverted, implying that the effective neutrino mass 
m l3/3 = 50 meV is well within the reach of next-generation neutrinoless double-beta decay experiments. 

Second, either the renormalization-group running effects of neutrino mixing parameters or the lepton 
flavor effects in leptogenesis has been ignored in the previous studies. Moreover, in order to stabilize the 
Higgs boson mass, the lightest heavy Majorana neutrino mass should be light enough M 1 < 4 x 10' GeV, 
which contradicts with the requirement M ] ~ 10 13 GeV for a successful leptogenesis. It is interesting to 
revisit this economical model by considering RG running effects, lepton flavor effects in leptogenesis and 
a non-hierarchical mass spectrum of heavy Majorana neutrinos. 

In this work, taking account of the RG running effects on neutrino mixing parameters, we have consoli¬ 
dated the conclusions reached in Ref. [13| and demonstrated that only four patterns B 1; B 2 , B 4 , and B 5 in 
Eq. (HD in the 10 case are allowed by current neutrino oscillation data. This generalization is important for 
the MSSM with a large value of tan /5, where the RG running effects are significant. It has been found that 
the determination of neutrino mass ordering and the observation of neutrinoless double-beta decays will 
provide critical evidences to verify or disprove these four patterns. Furthermore, the octant of 0 23 and the 
CP-violating phase 6 will be measured in future long-baseline neutrino oscillation experiments, and then 
can be used to further distinguish between B 4 (or B 4 ) and B 2 (or B 5 ). If the baryon number asymmetry 
is interpreted via leptogenesis mechanism, the relative sign of low-energy CP violation (i.e., the Jarlskog 
invariant J oc sin 5) to the high-energy CP violation (i.e., the CP asymmetry e 4 in N 1 decays) serves 
as a discriminator for B 4 (B 2 ) and B 4 (B 5 ). The most important formulas for four viable patterns are 
collected in Table 2. If the naturalness criterion is applied to the FGY model, only the nearly-degenerate 
mass spectrum of heavy Majorana neutrinos with a mass degeneracy of A ~ 1CP 7 is allowed, and resonant 
leptogenesis becomes responsible for the baryon number asymmetry. 

The FGY model actually exemplifies the idea of Occam’s razor, which cuts away unnecessary free 
parameters and renders the model to be most economical and testable. If one of four viable patterns of 
the flavor structure is singled out by future neutrino experiments, we should go further to identify the 
underlying symmetries and explore the true dynamics for neutrino masses and lepton flavor mixing. 
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Table 2: Collection of important formulas for the four viable patterns. 
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